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ABSTRACT 

With the discovery of Y dwarfs by the WISE mission, the population of field brown dwarfs 
now extends to objects with temperatures comparable to those of Solar System planets. To 
investigate the atmospheres of these newly identified brown dwarfs, we have conducted a pilot 
study monitoring an initial sample of three late T-dwarfs (T6.5, T8 and T8.5) and one Y-dwarf 
(YO) for infrared photometric variability at multiple epochs. With /-band imaging, each target 
was observed for a period of 1.0 h to 4.5 h per epoch, which covers a significant fraction of 
the expected rotational period. These measurements represent the first photometric monitoring 
for these targets. For three of the four targets (2M1047, Ross 458C and WISE0458), multi¬ 
epoch monitoring was performed, with the time span between epochs ranging from a few 
hours to ~2 years. During the first epoch, the T8.5 target WISE0458 exhibited variations 
with a remarkable min-to-max amplitude of 13 per cent, while the second epoch light curve 
taken ~2 years later did not note any variability to a 3 per cent upper limit. With an effective 
temperature of ~600 K, WISE0458 is the coldest variable brown dwarf published to-date, and 
combined with its high and variable amplitude makes it a fascinating target for detailed follow¬ 
up. The three remaining targets showed no significant variations, with a photometric precision 
between 0.8 and 20.0 per cent, depending on the target brightness. Combining the new results 
with previous multi-epoch observations of brown dwarfs with spectral types of T5 or later, 
the currently identified variables have locations on the colour-colour diagram better matched 
by theoretical models incorporating cloud opacities rather than cloud-free atmospheres. This 
preliminary result requires further study to determine if there is a definitive link between 
variability among late-T dwarfs and their location on the colour-colour diagram. 
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1 INTRODUCTION 

Ultracool dwarfs, spanning the L, T, and recently discovered Y 
dwarf (Kirkpatrick et al. 1999; Cushing et al. 2011) spectral types, 
provide a link between the coolest stars, giant planets in our So¬ 
lar System, and exoplanets. Without sufficient mass for nuclear 
fusion (e.g Hayashi & Nakano 1963), brown dwarfs cool mono- 
tonically over time, causing changes in the chemical and physical 
processes responsible for sculpting the emergent spectra of their 
atmospheres. The formation and dissipation of dusty condensate 
clouds are key components of theoretical models developed to ex¬ 
plain the fluxes and spectral features of brown dwarfs (e.g Allard 
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et al. 2001; Marley et al. 2002; Burrows, Sudarsky & Hubeny 2006; 
Helling & Woitke 2006). Early T dwarf atmospheric models pre¬ 
dicted that once the clouds from the L/T transition sink below the 
photosphere, the subsequent T-sequence should remain cloud-free 
(Marley et al. 2002). T-dwarfs however, appear to deviate from the 
expected cloud-free atmosphere colour as they cool (Figure 1; blue 
curve) and become progressively redder. This phenomenon can be 
best explained by the formation of sulfide and alkali salt clouds 
as the brown dwarf cools (Fodders & Fegley 2006; Visscher, Fod¬ 
ders & Fegley 2006; Morley et al. 2012). With a three-dimensional 
treatment of the atmospheric dynamics, recent models have sug¬ 
gested that large-scale temperature variations may be present near 
the photospheres of brown dwarfs and that these regions with dif¬ 
ferent temperatures may result in light curve variability (Showman 
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Figure 1. The L (red points), T (cyan points), and Y (purple points) brown 
dwarf sequence on an MKO J,H and iF2 {WISE 4.6p.m) colour-colour di¬ 
agram. The blue solid line represents a cloudless atmosphere showing the 
change in photometric colour with temperature (Saumon et al. 2012). The 
dashed red line and solid red line represent atmosphere models with an in¬ 
creasing cloud opacity from Morley et al. (2012). Magnitudes for the L, T & 
Y dwarfs are from Dupuy & Liu (2012); Kirkpatrick et al. (2012); Leggett 
et al. (2015). Details of the full sample are given in Table 1. 



& Kaspi 2013; Robinson & Marley 2014). Zhang & Showman 
(2014) also discuss how atmospheric turbulence and the resulting 
vortices and/or zonal jets can lead to photometric variability up to 
several percent over a wide range of timescales. 

Photometric monitoring of brown dwarfs provides a method 
to search for evidence and evolution of cloud features, storms, or 
activity that cause surface brightness differences (e.g. Koen et al. 
2005; Artigau et al. 2009; Buenzli et al. 2014; Radigan et al. 2014; 
Wilson, Rajan & Patience 2014). The strength, persistence, and 
wavelength dependence of the variations reveal the underlying at¬ 
mospheric processes of cloud formation, dissipation, and dynamics 
of the atmospheres (e.g. Showman & Kaspi 2013; Zhang & Show¬ 
man 2014). Over the past several years, we have conducted a large- 
scale photometric variability survey of brown dwarfs - the Brown 
dwarf Atmosphere Monitoring (BAM) program. The initial BAM 
study (Wilson, Rajan & Patience 2014) covered 69 targets span¬ 
ning the L0-T8 sequence, and detected multiple late-T variables. In 
this paper, we present the results of the second component of the 
BAM project that is designed to search for photometric variability 
over several epochs in four ultracool hrown dwarfs at or near the 
T/Y spectral boundary. The properties of the sample of T/Y dwarfs 
are summarised in Section 2. The near-infrared imaging observa¬ 
tions are described in Section 3. In Section 4, the data reduction 
and analysis required to construct the light curves is detailed. Re¬ 
sults for each target are given in Section 5, followed hy the discus¬ 
sion in Section 6 that includes a comparison with samples of higher 
temperature L and T brown dwarfs and theoretical models of cool 
atmospheres. 


2 THE BAM-II SAMPLE 

The sample for this pilot study consists of four ultracool field brown 
dwarfs which span the late-T to early-Y spectral types (see Figure 
1). Table 1 reports the target names, coordinates, /-hand magni¬ 
tudes, spectral types, effective temperatures and distances. The tar¬ 
gets were discovered as part of large-scale surveys with 2MASS 
(Burgasser et al. 1999), UKIDSS (Goldman et al. 2010), and WISE 
(Cushing et al. 2011). For this initial study, the targets include ob¬ 
jects with a range of properties, including youth (West et al. 2008), 
radio emission (Route & Wolszczan 2012; Williams, Berger & Za- 
uderer 2013), and binarity (Burgasser et al. 2012). 

The colour-colour diagram using near-infrared (/,//) and 
WISE (W2) filters in Figure 1 shows the variation in colours across 
the T and Y spectral types. Over-plotted is the clear atmosphere 
track from Saumon et al. (2012) and theoretical models of brown 
dwarf atmospheres that include the effects of emergent sulfide 
clouds in the photosphere from Morley et al. (2012). Clouds with 
different values of the parameter have different cloud proper¬ 
ties; low indicates smaller grain sizes and higher optical depth. 
Assuming that these clouds are patchy, we may expect to see more 
variability for the objects with colours consistent with more opti¬ 
cally thick (low /sed) clouds. Our targets span a range of colours 
and model tracks, ranging from completely cloud-free (blue line), 
to low /se 4 (solid red line) atmospheres with optically thick clouds, 
smaller grains, and larger vertical extent. 


2.1 2MASSW J1047539-I-212423 

The brown dwarf 2MASSW J1047539-I-212423 (2M1047) was 
identified with multi-epoch 2MASS images and classified as a 
T6.5 dwarf with Keck near-lR spectroscopy (Burgasser et al. 1999, 
2002). Using a parallax-based distance of 10.3 pc an effective 
temperature of ~900 K was inferred (Vrba et al. 2004). 2M1047 
is notable for measurement of radio variability, with bursts mea¬ 
sured at a frequency of 4.75 GHz (Route & Wolszczan 2012), and 
quasi-quiescent emission at 5.8 GHz (Williams, Berger & Zaud- 
erer 2013), making it the coolest brown dwarf with measured radio 
emission. Highly circularly polarised radio bursts lasting ~100 s 
were detected in three of the fifteen observations with Arecibo each 
with a cadence of 0.1 s over ~2 h. The radio observations were car¬ 
ried out over the course of one year, indicating a persistent source 
(Route & Wolszczan 2012). The quasi-quiescent emission was of 
longer duration (~40 min), but also two orders of magnitude fainter 
than the radio bursts. The target was monitored for variability in 
the J and H bands by Artigau, Nadeau & Doyon (2003) but did 
not exhibit any signs of variation. No contemporaneous optical or 
near-infrared photometric observations were recorded during either 
of the radio campaigns. 


2.2 Ross 458C 

Ross 458C is the 102 arcsec common proper motion substellar 
companion to the stellar binary Ross 458AB (Goldman et al. 2010), 
composed of two M-stars with a projected separation of ~5 au 
(Heintz 1994). The companion was detected in the UKIRT Infrared 
Deep Sky Survey (UKIDSS; Lawrence et al. 2007) and broad¬ 
band/methane filter photometry identified Ross 458 C as a late- 
T dwarf (Goldman et al. 2010). Subsequent near-infrared spec¬ 
troscopy determined a spectral type of T8 and effective tempera¬ 
ture of ~650 K (Burgasser et al. 2010). Similar assessments are 
reported in Burningham et al. (2011) and Cushing et al. (2011). 
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Table 1. Target list 


Target 

RA 

Dec 

Jmko 

(mag) 

SpT 

Temperature 

(K) 

Distance 

(pc) 

References 

2MASSW J1047539+212423 

10:47:53.8 

+21:24:23 

15.46 

T6.5 

~900 

10.6 

1,2,3 

Ross 458C 

13:00:41.7 

+12:21:14 

16.71 

T8 

~650 

11.7 

4, 5,6 

WISEP J045853.89+643452.9AB 

04:58:53.9 

+64:34:51.9 

17.47 

78.5“ 

-600“ 

14 

7,6 

WISEP J173835.53+273258.9 

17:38:35.52 

+27:32:58.9 

20.05 

YO 

430^50 

9.8 

11,9,6 


“ assigned value is for the binary and not the individual components. 

References: [1] Burgasser et al. (1999), [2] Burgasser et al. (2002), [3] Vrba et al. (2004), [4] Goldman et al. (2010), [5] Burgasser et al. (2010), [6] Dupuy & 
Kraus (2013), [7] Mainzeret al. (2011), [8] Delorme et al. (2010), [9] Cushing et al. (2011), [10] Liu et al. (2011) [11] Leggett et al. (2013). 


The distance to the Ross 458 system is 11.7 )*(q2q pc (Dupuy & 
Kraus 2013). The stellar pair in the system provides a means to es¬ 
timate the age through measurements of stellar activity. The age for 
a field brown dwarf is hard to constrain. Based on the strength of 
Ha emission, the level of variability (West et al. 2008) and space 
motion (Montes et al. 2001) of Ross 458AB, the age of the system 
has been estimated at <1 Gyr (e.g. Burgasser et al. 2010; Bum- 
ingham et al. 2011). Given the youth of the system, Ross 458C is 
predicted to be very low mass (5-20 Mjup; Burningham et al. 2011), 
which overlaps with planetary mass regime. Ross 458C is therefore 
a benchmark object for the investigation of both brown dwarf and 
exoplanet atmospheres. 

The atmosphere of Ross 458C measured by Burgasser et al. 
(2010) with near-IR spectroscopy reveals evidence of low surface 
gravity and the authors were able to better fit the near-infrared spec¬ 
trum with a cloudy atmosphere compared to cloudless. Burgasser 
et al. (2010) initially proposed the clouds in the atmosphere to 
be the reemergence of iron and silicate clouds, but more recently 
Morley et al. (2012) showed that models including salt and sulfide 
clouds fits the data better. Models incorporating both sulfide clouds 
and non-equilibrium chemistry have not yet been applied to the ob¬ 
servational data. 


2.3 WISEP J045853.89+643452.9AB 

WISEP J045853.89-h 643452.9AB (WISE0458) was the first ultra¬ 
cool brown dwarf discovered by the Wide-field Infrared Survey Ex¬ 
plorer satellite (WISE; Wright et al. 2010) in its search for the cold¬ 
est brown dwarfs in the solar neighborhood (Mainzer et al. 2011). 
Comparison of the medium resolution near-infrared spectrum of 
WISE0458 with a grid of cloudless models suggested a very cool 
effective temperature of ~600 K (Mainzer et al. 2011). High an¬ 
gular resolution imaging with the Keck laser guide star AO sys¬ 
tem revealed that the system has a binary companion at a separa¬ 
tion of ~0”5 and with a magnitude difference of ~1 mag (Gelino 
et al. 2011). Recent parallax measurements estimated the objects 
distance to be 14^j pc (Dupuy & Kraus 2013). AO-assisted spa¬ 
tially resolved spectroscopy confirmed both objects are very late T 
dwarfs near the T/Y boundary; from the resolved spectra, the pri¬ 
mary spectral type is T8.5 and the secondary spectral type is T9.5 
(Burgasser et al. 2012). We do not resolve the individual compo¬ 
nents in this study. 

2.4 WISEP J173835.53+273258.9 

The final target, WISEP J173835.53-I-273258.9 (WISEI738), is 
among the first Y dwarfs detected by the WISE satellite (Cush¬ 


ing et al. 2011) and is still only one of 18 known Y dwarfs (Kirk¬ 
patrick et al. 2012; Leggett et al. 2015). With a spectral classifi¬ 
cation of YO and an effective temperature of 430^4 q K (Dupuy & 
Kraus 2013), this object is one of the coldest brown dwarfs dis¬ 
covered. WISE1738 has been selected as the spectral standard for 
the YO class (Kirkpatrick et al. 2012). At these low temperatures 
water clouds are expected to form (Burrows, Sudarsky & Lunine 
2003), potentially leading to variability in the emergent flux. Inter¬ 
estingly, WISE1738 might show long period photometric variabil¬ 
ity. The object was originally measured to have a T-band magnitude 
of 19.51 ± 0.08 mag in Kirkpatrick et al. (2012), which was subse¬ 
quently measured to be 20.05 ± 0.09 mag in the recent Leggett et al. 
(2013) study. We adopt the latter magnitude for the purposes of this 
study. The Leggett et al. (2013) study does note the inconsistency 
in the original Kirkpatrick et al. (2012) photometry and suggests 
that the difference could be due to corrupted Palomar WIRC data. 
We adopt the Leggett et al. (2013) magnitude for the purposes of 
this study. 


3 OBSERVATIONS AND DATA REDUCTION 

This pilot study to monitor brown dwarfs at the T and Y dwarf 
boundary was initiated at the MMT observatory, with the first epoch 
of data for each of the four targets being obtained there. Follow-up 
observations were taken at three different observatories including 
the Canada-France-Hawaii Telescope (CFHT), UK Infrared Tele¬ 
scope (UKIRT), and the New Technology Telescope (NTT). The 
details for each of the targets is provided in the observing log pre¬ 
sented in Table 2. 

3.1 MMT 

Observations of the entire target dataset were taken with the SAO 
Widefield InfraRed Camera (SWIRC; Brown et al. 2008) at the 
6.5 m MMT Observatory in Arizona, on the 12th and 13th of March 
2012. The camera has an engineering grade 2048 x 2048 HAWAII- 
2 HgCdTe array, with a plate scale of 0.15 arcsec pixeU*, cor¬ 
responding to an on-sky field-of-view of 5.12 x 5.12 arcmin. We 
employed the T-band filter on SWIRC, which closely matches the 
Mauna Kea Observatory (MKO) T-band filter (Tokunaga & Vacca 
2005). This filter was selected since the largest amplitude varia¬ 
tions in known brown dwarf variables occur in the T-band (Artigau 
et al. 2009; Radigan et al. 2012). The observing strategy involved 
maintaining the target on a single pixel over an ~20 min timescale 
with a four-point dither pattern for the purpose of sky subtraction; 
this sequence was repeated over a ~1 - 4 hour time period, as sum¬ 
marized in Table 2. All of the raw images were calibrated using 
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Table 2. Observing log 


Target 

Telescope 

Date 

At 

(h) 

texp 

(s) 

Seeing" 

(arcsec) 


MMT 

2012-03-12T04:12:25 

1.0 

30 

0.6 

2M1047 

MMT 

2012-03-12T09:37:52 

1.9 

30 

0.8 


NTT 

2014-05-15T23:16:45 

2.8 

10 

0.8 


MMT 

2012-03-12T07:34:24 

1.5 

30, 60 

0.9 


UKIRT 

2014-04-23T07:02:54 

2.5 

10 

0.9 

Ross 458C 


NTT 

2014-05-16T23:42:38 

5.3 

20,40 

0.9 


NTT 

2014-05-20T04:18:56 

5.1 

40, 50, 70, 90 

0.9 


MMT 

2012-03-13T02:49:59 

3.7 

45 

1.0 

WTSE0458 


CFHT 

2014-03-21T05:50:09 

2.7 

60 

1.0 

WISE1738 

MMT 

2012-03-13T10:04:48 

2.3 

60 

0.8 


report the median seeing for each target. 


median combined darks and flat-field images. The exposure times 
ranged from 10s to 60s, depending on the target brightness, and the 
per-frame overhead for the detector was ~5 s. The SWIRC detector 
has a scattered light artifact in the lower-left quadrant of the detec¬ 
tor which was removed by using a high-pass filter. The SWIRC 
H2RG-detector is linear to within 0.1 per cent up to 40000 DN and 
the exposure times were set to ensure that the target brown dwarfs 
were maintained well below the levels approaching the non-linear 
regime of the detector. Additionally, comparison stars with peak 
fluxes greater than 40000 DN were rejected from the analysis. 

3.2 CFHT 

We obtained ~3 h of data on the T8.5 binary brown dwarf 
WISE0458 using the Wide-field InfraRed Camera (WIRCAM; 
Puget et al. 2004) on the 3.5 m Canada France Hawaii Telescope 
(CFHT). With a plate scale of 0.3 arcsec pixeP*, corresponding to 
an on-sky field-of-view of 20 x 20 arcmin, there were several tens 
of similar-brightness reference stars within the field of each target 
for the differential photometry calculation. The observations were 
carried out in a queue-based observing mode in T-band with a me¬ 
dian seeing of ~ 1" through most the sequence. The images were 
obtained in a staring mode, at a 60 s cadence. Per-frame overhead 
for the detector was ~6.5 s. The data was dark and flat field cal¬ 
ibrated using the T‘iwi automatic data pipeline from CFHT. The 
data was then sky subtracted in similar manner as described above, 
using a median filter (where the stars were first masked) to generate 
individual sky frames. The CFHT detector is linear to better than 
a per cent up to 5000 DN. The target and all the reference stars 
were exposed to less than this limit. Comparison stars brighter than 
5000 DN were rejected in the analysis. 

3.3 UKIRT 

The second epoch on Ross 45 8C was obtained on 23 April 2014 us¬ 
ing the Wide-Field infrared Camera (WFCAM; Casali et al. 2007), 
an infrared wide-field camera on the 3.8 m UK Infrared Telescope 
(UKIRT). WFCAM has four 2048x2048 array detectors with a 
pixel scale of 0.4 arcsec pixeU*, corresponding to an on-sky field- 
of-view of 0.21 X 0.21 degrees. The observations were carried out 
with the MKO /-band filter, with seeing of ~ 0.9 arcsec during the 
sequence. A five-point dither pattern was used, with each individ¬ 
ual image having an exposure time of 10 s. Per-frame overhead for 
the detector was ~1.3 s. The data were calibrated, and the differ¬ 
ent dithers were combined using the a dedicated pipeline developed 


by the Cambridge Astronomy Survey Unit (Irwin 2008). The WF¬ 
CAM detectors are linear to better than a per cent within 40000 DN 
and care was taken to ensure that neither the target nor the reference 
stars used in the reduction were approaching this limit. 

3.4 NTT 

2M1047 and Ross 458C were observed between the 15th and 23rd 
of May 2014 with Son of ISAAC (SofI; Moorwood, Cuby & Lid- 
man 1998) mounted on the 3.6 m NTT (New Technology Tele¬ 
scope). The observations utilized the wide-field imaging mode of 
SofI, with a plate scale of 0.28 arcsec pixeU*, corresponding to an 
on-sky field-of-view of 4.92 x 4.92 arcmin. Per-frame overhead for 
the detector was ~7.5 s. The NTT detector is linear to > 1.5 per cent 
when objects have less than 10000 DN peak flux. The observations 
were all obtained using the filter, with a two-point AB-AB nod¬ 
ding pattern based on recommendations from the instrument scien¬ 
tist. The flux of the target, and reference stars within the field, was 
kept below 10,000 DN to minimise the effects from the detector 
non-linearity. To limit systematics in the data we used a two-point 
nod which permitted an accurate estimate of the sky background for 
the targets. For each object, data reduction consisting of correcting 
for the dark current and division by a flat field and sky subtraction 
were applied. 


4 LIGHT CURVE GENERATION AND IDENTIEICATION 
OE VARIABILITY 

We performed aperture photometry on the calibrated and aligned 
images from each of the observatories, using the APPHOT pack¬ 
age in IRAF^. For each of the targets in our study we computed 
the aperture photometry for a range of aperture sizes ranging from 
radii of 0.6 - 2.0 times the full width half maximum (FWHM), 
and found that an aperture of 1.0 x FWHM provided the highest 
signal-to-noise ratio (SNR) across the full sample. Data taken dur¬ 
ing periods of poor conditions were not included in the analysis; 
the frame selection caused the gap between epoch 1 and 2 from 
2M1047 and limited the W1738 observation period. A catalogue 
of all the stars in the field of view was generated, which typically 
included 10 ~ 50 stars, from which the 15 reference stars most sim¬ 
ilar in brightness to the target were selected. We chose the similar 
brightness stars over the brightest stars for two reasons: to limit 
non-linearity effects, and to ensure that non-astrophysical varia¬ 
tions due to weather were more accurately duplicated by the ref¬ 
erences. The limited field-of-view of the different detectors meant 
that there were not any objects of the same late-T spectral types to 
be used as references. 

The final target light curves were generated by the photo¬ 
metric pipeline developed as part of the BAM-I survey (Wilson, 
Rajan & Patience 2014), modified to measure the photometry on 
each individual image. Measuring the photometry on individual 
images and taking the median value or forming a median image 
and measuring photometry on the median combined image did not 
result in significantly different light curves. The determination of 
whether the observed flux variation in the target light curve was 

* IRAF is distributed by the National Optical Astronomy Observatories, 
which are operated by the Association of Universities for Research in 
Astronomy, Inc., under cooperative agreement with the National Science 
Foundation. 
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Figure 2. Detrended multi-epoch light curves for each of the targets with the master reference light curve (light blue). Top-left three panels are the different 
epochs for 2M1047 (black), with the date of observation and telescope noted in the panel. Bottom-left two panels are the same for W0458 (red). Top-right 
four panels are the multi-epoch observations for Ross 458C (orange). The bottom-right panel is the single epoch obtained for the YO brown dwarf WISE1738 
(green). 


due to an astrophysical process, required that the objects have a p- 
value 5 per cent. The p-value is defined as the probability that 
the final target light curve is the same (p-value > 10%) or differ¬ 
ent (p-value ^ 5%) from the master reference light curve i.e. the 
median combination of all the reference light curves. Calculation 
of the p-value statistics for each of the light curves was carried out 
on weighted-mean combined data points where the errors are the 
1 C7 scatter in each bin, rather than on the light curves composed of 
the individual photometry points. This was done to ensure that the 
statistics were not biased by outlier data points with small errors. 


5 RESULTS 

The T-band light curves for each target generated using the obser¬ 
vations at each epoch are shown in Figure 2, plotted alongside the 
master reference light curve used to determine whether variations 
were observed. The master reference light curve is generated by 
median combining all the reference star light curves and shows any 
residual trends common to all light curves. A summary of the re¬ 
sults of the light curve analysis, including the peak-to-peak ampli¬ 
tude of any variation and the associated p-value, is given in Table 3. 
Each light curve is normalised, and variability is only investigated 
over the time scale of an individual epoch; no attempt to distinguish 
photometry variations between the individual epochs was made. 


5.1 2M1047 

The time-series photometric measurements for the T6.5 dwarf 
2M1047 are shown in the three top-left panels of Figure 2. The 
first two epochs are separated by ~3.5 h and the final epoch of 
data was taken 2.2 yr later. None of the three light curves for this 
target show statistically significant variability, and these results are 
similar to a previous y-band monitoring with the 1.6 m telescope 
of the L’observatoire du Mont-Megantic by Artigau, Nadeau & 
Doyon (2003). Although there were no contemporaneous radio ob¬ 
servations at the time of our y-band imaging, previous measure¬ 
ments of 2M1047 has shown variability at radio frequencies (e.g. 
Route & Wolszczan 2012). The radio bursts occurred over a time 
period of ~100 s which is less than a single binned data point in 
the light curves. Radio bursts were recorded three times out of the 
fifteen observations (2 h each) made with Arecibo, but the current 
near-infrared data show no intensive brightening over comparable 
13 h time intervals. In addition to the intense radio bursts, quasi- 
quiescent fluctuations in the radio emission with a timescale of 
'-^45 min have been reported (Williams, Berger & Zauderer 2013), 
however these variations were a factor of a hundred fainter than the 
large bursts described in Route & Wolszczan (2012). Contempo¬ 
raneous radio and near-IR observations are required to search for 
any correlation between radio bursts and variations in photospheric 
flux. 
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Table 3. Summary of BAM-II variability study 


Object 

Spectral Type 

Telescope 

Duration 

Binned Images* 

References 

V 

xl 

p-value 

Amplitude/Limit** 




(h) 





(%) 

(%) 



MMT 

0.95 

12 

7 

6 

3 

20 

<0.8 

2M1047 

T6.5 

MMT 

1.88 

19 

7 

10 

1.92 

10 

< 1.1 



NTT 

2.75 

10 

5 

15 

0.79 

98 

< 1.0 



MMT 

1.78 

7 

6 

11 

0.58 

10 

< 1.1 

Ross 458C 

T8 

UKIRT 

2.53 

5 

10 

7 

8 

0.36 

10 

94 

<0.8 

NTT 

5.29 

9 

13 

0.41 

< 1.8 



NTT 

5.11 

10 

7 

13 

1.12 

93 

<2.1 

WISE0458 

T8.5 

MMT 

CFHT 

3.68 

11 

8 

14 

12 

11.37 

0 

84 

13±3 

2.74 
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9 

0.21 

<2.6 

WISE1738 

YO 

MMT 

2.39 
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10 

5 

0.9 

15 

<20.3 


Notes: *The average number of individual images used to generate the weighted-mean light curves. ** For the constant light curves the limit is the l-ff 
photometric error of the target light curve. For the variable epoch of W0458, the amplitude is the peak-to-trough value of the best fit sine curve to the data. 


5.2 Ross 458C 

The four light curves for Ross 45 8C are given in the top-right pan¬ 
els of Figure 2. The time spans between observations ranged from 
days to yr, with 2.1 yr from the first to second epoch, 23.7 days be¬ 
tween the second and third epoch, and 3.2 days between the third 
and fourth epoch. No statistically significant variations were de¬ 
tected at any of the four epochs, and the limits on detectable am¬ 
plitudes ranged from 0.8 to 2.1 per cent. The results suggest that 
there are no large and persistent storm features that would induce 
rotationally modulated brightness changes, or that the system is 
viewed pole-on. The Ross 458C data form the most comprehen¬ 
sive monitoring of the atmosphere of a brown dwarf that serves as 
an exoplanet analogue. Ross 458C represents the later stage of at¬ 
mosphere evolution compared to the younger imaged exoplanets 
such as HR8799 d and j3 Pic b that may have similar masses but 
are substantially warmer due to their younger ages. Ross 458C also 
occupies an intermediate location in the colour-magnitude diagram 
between the youngest directly imaged planets and the older, cooler 
GJ504 b exoplanet (Kuzuhara et al. 2013), but is far less technically 
challenging to monitor because of the wider angular separation be¬ 
tween Ross 458AB and Ross 458C. 



Time (hours) 


Figure 3. WISE0458 light curve (red circles) for the variable MMT epoch 
obtained on the 13th March 2012. Also plotted is the best fit sine wave (red 
dashed line) to the W0458 light curve resulting in a period of ~3 h with a 
variable min-to-max amplitude of 13 per cent. The master reference light 
curve is plotted with the light blue points. 


5.3 WISE0458 

The bottom left two light curves in Figure 2 show the dramatic 
difference between the two epochs of observations for the binary 
brown dwarf WISE0458. In the first epoch, the target is highly vari¬ 
able compared to the master reference, with a min-to-max ampli¬ 
tude of ~ 17 per cent. The measured variability of WISE0458 is the 
second highest amplitude brown dwarf observed to date in a brown 
dwarf, with only 2M2139 exhibiting greater variability (Radigan 
et al. 2012). The first epoch of WISE0458 data exhibits a periodic 
pattern and we fit a series of pure sine waves of different amplitudes 
and periods to the light curve. The best fit periodic signal with an 
amplitude of 13.2 per cent and a period of 3 h is shown in Figure 3. 
The second epoch data on WISE0458 taken after a gap of ~2 yr 
lacks any detectable variations. Further monitoring of WISE0458 
will determine if the large amplitude variations recur. 


5.4 WISE1738 

The light curve for the YO target WISE1738, shown on the bottom 
right of Figure 2, is impacted by the large uncertainties and the lim¬ 
ited number of data points. During initial planning of the observing 
run at the MMT, WISE1738 was not expected to be the faintest tar¬ 
get in our sample, however, as noted in Section 2, the photometry 
of this target was re-estimated to be 0.5 magnitudes fainter than 
originally measured (Leggett et al. 2013). This intrinsic faintness 
combined with variable seeing during the observation resulted in 
low signal-to-noise ratio detections of the target and, consequently, 
an extremely noisy light curve. The light curve does not show any 
statistically significant variations, but the limitations in the photo¬ 
metric precision restrict the interpretation of the source as a con¬ 
stant with a large limit on possible variations of < 20 per cent. 
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Figure 5. Plot summarizing literature multi-epoch variables and constants. Only objects with multiple epochs are plotted in the figure, with targets coming 
from (Artigau et al. 2009; Buenzli et al. 2012; Clarke et al. 2008; Girardin, Artigau & Doyon 2013; Khandrika et al. 2013; Koen 2013; Metchev et al. 2013; 
Radigan et al. 2014; Wilson, Rajan & Patience 2014). The BAM-II targets are indicated with black symbols for 2M1047, orange symbols for Ross 458C, 
red symbols for WISE0458 and green symbols for WISE1738. Each vertical line corresponds to a unique object, and solid lines indicate brown dwarfs that 
remained consistently variable or constant at more than one epoch, while the dashed lines identify the objects that switched between variable and constant 
states. The shaded regions indicate the L/T transition (pink) and TA^ boundary (gray) regions. 



Spectral Type 

Figure 4. Spot fractional coverage required for a 2 per cent amplitude J- 
band light curve, plotted as a function of spectral type. The fractional cover¬ 
age was estimated using the solar metallicity, log(g) = 5.0 synthetic spec¬ 
tra within the BT-Settl model grid (Allard, Homeier & Freytag 2011). 
For each Tgff within the grid, the fractional projected surface area at a 
higher/lower temperature (ATeff = ±300) required to increase/decrease the 
7-band flux by 2 per cent was estimated. This first-order approximation is 
analogous to the variability induced by a region of higher or lower temper¬ 
ature rotating in and out of view. 

6 DISCUSSION 

The BAM-II study was designed as a pilot program to investigate 
the coolest brown dwarf atmospheres using multi-epoch photomet- 



J-H (mag) 

Figure 6. Colour-colour diagram showing the sample (large filled stai's) and 
the Y dwarfs (purple) along with variable (red) and constant (blue) brown 
dwarfs with spectral type greater than T5 from Figure 5. Figure 1 shows 
the full sample. The blue line is theoretical track for clear atmosphere from 
Saumon et al. (2012) and the red lines indicating cloudy atmospheres are 
from Morley et al. (2012). All of the models assume a surface gravity of 
log(g) = 5. 
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ric monitoring as a probe of the dynamics and surface brightness 
variations. To place the BAM-II results in the broader context of 
the brown dwarf population. Figure 5 combines the targets in both 
this study, and the previous BAM-I survey spanning the full L0-T8 
sequence (Wilson, Rajan & Patience 2014). Single epoch BAM-I 
variables are also included in the comparison. Additional epochs 
beyond the BAM-I measurements were provided by measurements 
from a recent large survey of L3-T9 brown dwarfs (Radigan et al. 
2014) and from the more focused studies compiled in the BAM-I 
paper (Artigau et al. 2009; Buenzli et al. 2012; Clarke et al. 2008; 
Girardin, Artigau & Doyon 2013; Khandrika et al. 2013; Koen 
2013; Metchev et al. 2013; Radigan et al. 2012). 

Nearly all of the amplitudes or limits shown in Figure 5 were 
measured in the T-band, however five objects include results from 
different filters, since this was the only way to include more than 
one epoch for those targets. Both the significant number of objects 
that switch between variable and constant states and the substantial 
range of amplitudes for the multi-epoch variables highlight the dy¬ 
namic and evolving nature of substellar atmospheres, as well as the 
need for multi-epoch monitoring. Currently, the best studied case 
is that of SIMP0136 with monitoring from 2008 to 2012 revealing 
a remarkable evolution in both the amplitude (from > 5 per cent to 
undetectable) and the shape (from sinusoidal to multi-component) 
of the light curve measured over several hours per epoch (Metchev 
et al. 2013). The BAM-I and BAM-II surveys have identified a set 
of targets spanning the full L-T sequence that warrant further mon¬ 
itoring. 

The variable and constant brown dwarfs with spectral types of 
T5 and later that are included in Figure 5 are plotted on an (H — W2) 
vs. (J — H) colour-colour plot in Figure 6 to compare with a set 
of theoretical models (Saumon et al. 2012, Morley et al. 2012). A 
lack of WISE photometry for some targets in Figure 5 limits the 
number of brown dwarfs from this study that can be plotted in Fig¬ 
ure 6. The distribution of the late-T and Y dwarf population in this 
colour-colour diagram (larger sample in Figure 1) cannot be recon¬ 
ciled with models of cloud-free atmospheres (see Figure 6), and a 
proposed explanation for the dispersion involves varying amounts 
of opacity parametrised in the value of /sed (Morley et al. 2012). 
The redder J — H colour of WISE0458 places the target close to 
the /sed = 4 model with intermediate level of cloudiness as indi¬ 
cated in the colour-colour diagram of Figure 6. The constant targets 
(blue circles) amongst the late-T objects with multi-epoch measure¬ 
ments presented in Figure 6, appear concentrated near the cloud- 
free model, and the variables (red circles) appear to have redder 
and thus cloudier atmospheres. Changes in the /sed values may be 
linked not only to the colour, hut also to the presence or absence 
of variability. Determining whether or not there is a link between 
the variability and location in the colour-colour diagram requires a 
larger set of late-T monitoring observations. 

Recent work on synthetic atmosphere models extending to the 
temperatures of late-T and Y dwarfs (900 — 400 K) has indicated 
that these brown dwarfs may have sulfide and alkali salt clouds in 
their photospheres (Morley et al. 2012). Such clouds peak in op¬ 
tical depth for objects with Teff ~ 600 K (T9) but persist in ob¬ 
jects from 900 K to under 400 K. Three-dimensional models that 
include radiative transfer and cloud formation have not yet been 
developed for these cool objects. Idealised 3D-circulation simu¬ 
lations by Showman & Kaspi (2013) suggests that brown dwarfs 
may have complex circulation patterns on regional and global 
scales, generated by the interaction of convective layers with the 
overlying stably stratified radiative atmosphere, potentially lead¬ 
ing to patchy cloud structure. As stated in the introduction, atmo¬ 


spheric circulation models including atmospheric turbulence pre¬ 
dict a range of variability amplitudes over multiple epochs for the 
emitted flux Zhang & Showman (2014). If cool brown dwarfs do 
have weather patterns causing heterogeneous cloud cover, the Mor¬ 
ley et al. (2012) models predict that they would show photometric 
variability in the near-infrared, predominantly in Y and /-bands. 
Based on the initial results from this study, there is an indication 
that the variables have redder colours and are concentrated in the 
region of the colour-colour diagram associated with the atmosphere 
models more likely to have weather patterns. This suggestive link 
between theory and observations will be investigated further with 
larger surveys for multi-epoch variability among the coolest brown 
dwarfs. 

Motivated by recent studies by Robinson & Marley (2014); 
Morley et al. (2014), we investigated what variability of different 
per cent amplitudes mean in terms of actual hot/cold spot fractional 
coverage on the brown dwarf photosphere. We ran a simple simula¬ 
tion using the BT-Settl model grid from Allard, Homeier & Frey- 
tag (2011), similar to what was done in Rostov & Apai (2013) but 
for a larger temperature range. In the simulation we estimated the 
spot coverage required to produce a particular amplitude of vari¬ 
ability as a function of the spectral type, where the spectral types 
are defined by their temperature. Figure 4 shows the results of a 
simulated 2 per cent variable, for which the patches have a ATeff 
of ±300K. The figure indicates that for different spectral types, the 
spot fractional coverage required to produce the same amplitude of 
variation can differ by several percent. And amplitude cutoffs to in¬ 
dicate “strong” or “weak” variations might not necessarily indicate 
higher and lower spot fractional coverage. 
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